Ferroelectricity due to orbital ordering in E-type undoped rare-earth manganites. 
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Aiming at understanding the origin of the electronic contribution to ferroelectric polarization in 
undoped manganites, we evaluate the Berry phase of orbital-polarizable Bloch electrons as an orbital 
ordering (00) establishes in the background of an antiferromagnetic E-type configuration. The onset 
of 00 is tuned by the Jahn- Teller (JT) interaction in a tight-binding model for interacting electrons 
moving along zigzag chains. A finite polarization is found as soon as the JT coupling is strong 
enough to induce 00, supporting the large electronic contribution predicted from first principles. 
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PACS numbers: 75.47.Lx, 77.80.-e, 75.25.Dk, 71.70.Ej 



The coexistence of ferroelectrically polarized state with 
long-range magnetic order is usually referred to as multi- 
ferroic behavior [lj]. Among other multiferroic materials, 
orthorhombic rare-earth manganites i?Mn03 (R belong- 
ing to lanthanum series) represent an important class of 
"improper multiferroics" Q , where the ferroelectric po- 
larized state is not only coexisting with, but also intrin- 
sically related to some kind of magnetic order. Several 
degrees of freedom (spin, charge, orbital, lattice) are ac- 
tive in these systems and are responsible for their very 
rich phenomenology. Starting from R =La and moving 
along the series to smaller ionic radius r^, the ground 
state changes from antiferromagnetic A-type (AFM-A) 
to E-type (AFM-E), through an incommensurate spi- 
ral structure for R =Tb, Dy[3j. The incommensurate 
magnetic phase of TbMnC>3 and DyMnC>3 is responsi- 
ble for the observed weak ferroelectric polarization^, 
which is ultimately driven by spin-orbit coupling via 
the Dzyaloshinskii-Moriya interaction 0, |(|. The rela- 
tivistic origin of ferroelectricity in these compounds is 
reflected in the weak P (smaller than 0.1/iC/cm 2 [|j]). 
On the other hand, also collinear AFM-E magnetic or- 
der has been predicted to induce a ferroelectric polariza- 
tion via an exchange-striction mechanism[7[, where the 
double-exchange interaction between Mn e g electrons in 
the symmetry-broken ft44- spin configuration is respon- 
sible for polar atomic displacements of oxygens bridg- 
ing Mn ions. In a following paper Q, this picture has 
been partially confirmed by means of density-functional 
theory (DFT) calculations, reporting P ~ Q^,C/cm 2 
for o— H0M11O3. The DFT analysis also pointed out 
that a comparable contribution to the total polariza- 
tion has a purely quantum mechanical origin, possibly 
related to the AFM-E-induced asymmetric electron hop- 
ping of orbitally polarized e x g states. The microscopic 
origin of such electronic contribution to P has been pro- 
vided in terms of maximally localized Wannier functions 
(WFs)@. The DFT study showed that WF centers in all 
AFM-E o - i?Mn0 3 (R =Ho, Er, Tm, Lu), are largely 
displaced with respect to corresponding ionic positions; 
thus, the total P arises from the sum of an ionic con- 
tribution (driven by the exchange-striction mechanism) 
and a purely electronic one@. 

From the experimental point of view, only a small 



polarization (P ~ 0.01/iC/cm 2 ) has been reported 
for o— H0M11O3 10], at odds with theoretical predic- 
tion. Only recently experimental results for AFM-E 
o— TmMn03 suggested that P could actually pass the 
1/j.C/cm 2 threshold [ll|. By the way, we notice that 
the thermodynamically stable phase for J?MnOa with 
rfl smaller than that of Dy is hexagonal rather than 
orthorhombic [l^: the desired perovskite structure can be 
obtained, e.g., by high-pressure synthesis, leading how- 
ever to polycrystalline samples that can limit the accu- 
racy in measuring the predicted properties. 

Aim of this Letter is to set a clear correspondence 
between the electronic polarization and the onset of or- 
bital ordering (00) on the background of the magnetic 
AFM-E configuration. The proposed picture is com- 
pletely new in the field of improper multiferroicity and 
relies on the orbital polarizability of electrons on a spe- 
cific magnetic background. We will introduce a simple 
model which is expected to reproduce the general features 
of the ground state of AFM-E orbital-ordered i?Mn03. 
The AFM-E configuration of ti g spins is treated as a 
constraint on the e g electrons motion, a reasonable as- 
sumption in the limit of very large Hund coupling. Its 
stability has been already discussed, e.g. in Refs. [!],[§], 
where the competition between kinetic energy (double 
exchange) and superexchange interactions between the 
Mn t2 g spins is discussed as a function of the ionic ra- 
dius rn (as pointed out recently, the Jahn- Teller (JT) 
interaction also may play a relevant role[13]). In limit 
of infinite Hund coupling, e g electrons can only hop be- 
tween sites with ferromagnetically aligned core-spins; this 
implies that hopping processes in the AFM-E phase are 
allowed only within one-dimensional zigzag chains of par- 
allel core spins [FigQJa)]. 00 may be stabilized by 
an electron-lattice JT interaction, which causes also the 
Bloch electrons within the chains to acquire a geomet- 
ric Berry phase arising from a conical intersection of the 
potential energy surfaces [l4j: we stress the fact, however, 
that also a correlation-mediated mechanism may stabilize 
OOflij. I n the framework of the modern theory of polar- 
ization, the appearance of this Berry phase can give rise 
to a polarization whose origin is purely quantum mechan- 
ical, being intimately connected to a current flow inside 
the bulk[16]. We will show that this is indeed the case, 
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but that the geometric phase of the 00 state is actually 
not enough, being the interplay with the specific topol- 
ogy induced by the underlying magnetic background the 
boost for ferroelectric polarization. 

Let us introduce the Hamiltonian which describes the 
motion of interacting e g electrons within zigzag chains: 



H = ~ E ^7' a \j a Yl + ^ E Ua 3 n « 
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+ E,,t [ 2 ( l2jT x j + q 3 jT ZJ ) + q\ + qlj] 
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The first term describes the electron transfer between 
nearest-neighbor Mn sites, where ajj, ■ creates a particle 
at site j in electronic states stemming from Mn +3 or- 
bitals d x 2_ y 2 (a) and d^ z 2 _ r 2 (/?) . The hopping ampli- 
tudes are t aa = 3io/4, tpp = to/4 and t a p = ±^3^/4 
(l7j . where t — (pda) 2 is the energy unit and the sign 
appearing in t a p depends on the hopping direction along 
the zigzag chain (i.e. t^g — —t v aB ). As we will see, 
this implies that at each site electrons gain a phase that 
depends on the orbital through which they pass. The 
second term describes the interorbital interaction, the 
only Coulomb interaction that is left between e g elec- 
trons in the infinite Hund coupling limit, while the last 
term accounts for the JT interaction with dimensionless 
Qisj = {k/g)Qvj, where Q V j are the Jahn- Teller-active 
modes[18j, k,g are respectively the elastic JT stiffness 
and the bare JT coupling, and Ejt — g 2 /{2k) is the 
static JT energy. r w = V 



- y^^'Vj is the or " 



bital pseudospin, being cr , the Pauli matrices, while 



n ii = a \j a ij are or bital density operators. We treat the 
on-site correlation in a mean-field framework, by lineariz- 
ing it in the optimum local basis in order to keep track 
of the orbital degrees of freedom[19], finding n aj n^j w 
(2nj (rij)- 2t zj (t zj ) - 2t xj (t xj ) - {n 3 ) 2 + (t zj ) 2 + (r xj ) 2 ) /4. 

As the phase change in hopping amplitudes plays a 
relevant role in establishing a ferroelectric polarization, 
we transform the e g electron basis through the unitary 
transformation Cj = (a a j + iapj)/V2 and dj = (a a j — 
i a/jjVv^tlJ]. In the new basis the Hamiltonian (TTJ) is 
rewritten as 



H = - 
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+ E JT ( q % + 4) - V - {{n 3 ) 2 - (r ZJ ) 2 - (r x ,) 2 ) 



with t = to/2 and s = e'^-Jl to/2, the phase 4>\%-j\ 
depending on the hopping direction as 4> y — — <fi x = ir/3. 
Therefore, the e g electrons pick up a phase change as 
they move between different neighboring orbitals. The 
local interaction Vj — \Vj \ e 1 ^ acquires a phase too flij . 
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where we set u = U/2Ejt, while (t^j) are the aver- 
aged pseudospin operators in the original basis to be 
self-consistently determined. As in undoped manganitcs 
there is only one e g electron per site, according to Ref. 
[H[ we can take qsj = q$ and q2j = {—l) ] q2 [we numeri- 
cally checked this assumption by minimization of Hamil- 
tonian ((2|) with respect to q a i\. The absolute value of the 
interaction potential is then found to be constant, while 
its phase changes within the chain with a periodicity 
equal to 7r. The Hamiltonian in momentum space can be 
rewritten then in the electron basis (c£, d\, <4 +7r , d\ +7r ), 
being k defined, with the unit cell chosen as shown in 
Fig[lla), in the reduced Brillouin zone —ir/2 < k < ir/2: 
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The hopping phase change enters in k± = k ± 7r/3, 
whereas the Fourier components of the local potential 
are Vo = |y|cos£_e^+ and V v = \V\i sin£_ e^+, be- 
ing ^± = (£i ± £,2)/2. We note that the system is al- 
ways in a band-insulating state, even when the JT cou- 
pling and U are set to zero. The band insulator is sta- 
bilized by the phase difference between the interorbital 
hopping amplitudes, as it happen s in the AFM-CE phase 
of half-doped manganites[l9|, [20|, and it is very robust 
due to its topological origin. However, the 00 emerges 
only when the JT and Coulomb interactions are suffi- 
ciently strong to induce a nonzero \V\. In order to in- 
spect the 00 in real space, we can evaluate the local 
average orbital occupancy p-yj = ((/^ ■</> ■), where 7J - = 
— sm(9jj/2) a a j + cos(6* 7J /2) apj and # T ; = 27r/3(47r/3) 



for 7 = 3x — r 2 (3y 



! ). In FiglUc) we show the 



evolution of p 3x 2_ r 2 on site 1 (equivalent to p 3y 2_ r 2 on 
site 2) as a function of JT interaction at U = and the 
00 pattern of alternating 3x 2 — r 2 /3y 2 — r 2 within the 
chain. We note that a similar result is obtained for the 
linear chain, where the hopping direction and the related 
phase 4>\i-j\ are unchanged, while the OO-induced phase 
change is retained [l4j]: in this case the band insulator has 
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FIG. 1: (a) Unit cell for the zigzag ferromagnetic chain in 
the AFM-E phase: the hopping direction changes periodi- 
cally within the chain as {..x,y,x,y..}. (b) Orbital densities 
Ps x 2_ r 2 , Ps y 2_ r 2 as a function of position within a chain at 
Ejt = 1, when OO is almost fully developed, (c) Average 
orbital density for d 3:E 2_ r 2-state at site 1 as a function of JT 
interaction (inset: total energy smooth evolution). Dotted 
line is the same quantity evaluated in a linear chain, when 
<P\i—j\ — const, (d) Same as (c), but for different values of 
interorbital Coulomb interaction. 



no topological origin and emerges in an orbitally ordered 
pattern because of the local interaction, being metallic 
for \V\ — [dotted line in FigQTc)]. On the other hand, 
the inclusion of U, albeit at a mean-field level, makes the 
OO onset more robust [FkjQT'd)]. 

We consider now the possible ferroelectricity. Follow- 
ing the prescription described in Ref. [2lj, we evaluate 
the polarization as 



P = — — lim Im In det S, 

Z7T L— >oo 



(5) 



where L is the chain length (periodic boundary condi- 
tions are assumed) and S is the overlap matrix, defined 
as S m<m >(k,k') = (^ m ,k\e~ l2 ^ x \^ m ',k')- The eigenvectors 
ipm k, where to is a band index, can be expressed as 



(U k (m,l)cl + U k (m,2)dl+ 
U k (m,3)cl +n + U k (m,4)4 + ) |0). (6) 



IV'fc.m) = 



Here Uk is the unitary matrix which diagonalizes the 
Hamiltonian Q at each k, and |0) is the vacuum state. 
By defining the position operator as x = y\ j (cj + 
dj dj), one sees that S elements vanish except when each 

pair of vectors k, k differs by an amount e — 2ir/L; then 
the determinant can be factorized into L small determi- 
nants whose dimension is equal to the number of occu- 
pied bands [2l|, giving det S = Uk det S(k, k + e), with 



the small overlap matrix being 

S m , m >(k,k + e) = U f k (l,m)U k+E (m' n ). (7) 

7=1,4 

In Figj2] we show the polarization evaluated through 
formula ([5]) for the ground state of our model: it remains 
equal to zero in the insulating phase as far as no OO 
is induced in the zigzag chain by the local interaction, 
then it rapidly increases, closely following the evolution 
of the occupied orbital density. The sign of P can be 
changed by rotating the OO pattern, in such a way that 
orbital state 3y 2 — r 2 (3x 2 — r 2 ) is occupied at site 1(2) 
instead of 3a; 2 — r 2 (3y 2 — r 2 ). It is worthwhile to notice 
that P is always zero in the orbital-ordered insulating 
phase found in the linear chain, even though the interac- 
tion phase £j has been shown to induce a Berry phase in 
the electronic Bloch function for a given k\l4- Indeed, 
this geometric phase takes into account the difference 
in the electron motion, clockwise or counterclockwise, 
around each site displaying OO within the linear chain; 
however, this difference sums up to zero when evaluated 
along the whole linear chain. Similarly, when the local 
interaction is not strong enough to induce OO along the 
zigzag chains, the difference in the electron motion due 
to direction-dependent interorbital hopping amplitudes, 
which induces a phase change in the electronic Bloch 
functions, does not give rise alone to any nonvanishing 
P. On the other hand, our result suggests that the inter- 
play between the phase changes induced by the OO and 
direction- dependent hopping amplitudes is responsible for 
the onset of a ferroelectric state, which has ultimately a 
topological origin. We stress that this polarization has a 
purely electronic origin, since the position of the ions is 
fixed within the chains, supporting the DFT calculations 
which report a large electronic contribution to P in the 
whole class of AFM-E o - i?Mn0 3 [l|. We can estimate 
the magnitude of the calculated electronic P for a realis- 
tic three-dimensional system multiplying it by a scale fac- 
tor ao/Vo, where ao is the Mn-Mn distance (correspond- 
ing to the lattice constant in our simplified model) and Vq 
is the unit cell volume. We estimate ao = 3.89 — 4 A and 
V = 244 - 220 A 3 for i?Mn0 3 [I|, implying P el of the 
order of ~ 10 \iC '/cm 2 . On the other hand, both Ejt /to 
and U/to are of the order of unity in manganites, being 
to = 0.1-0.5 eV, E JT ~ 0.25 eV and U > 1 eV^, sug- 
gesting that the OO found in our simple model is easily 
realized in undoped manganites. 

To get more physical insight, we evaluated the position 
of the WF centers in the chain. For this purpose, the 
construction of maximally localized WFs is not needed 
in one-dimensional systems. Indeed, they can be ob- 
tained as the eigenvalues of a matrix A, constructed as 
the product of the unitary parts of the S matrices along 
the fc-point string (by "unitary part" we mean the matrix 
product VW^ taken from the singular value decomposi- 
tion 5 = VXW^, where V and W are unitary and S is a 
diagonal matrix with nonnegative diagonal elements) [23| . 
We find that i) in the insulating phase with no OO, the 
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FIG. 2: Polarization in the one-dimensional zigzag chain as a 
function of Ejt; the electron charge e is taken equal to unity. 
Solid (dashed) line is evaluated for the 3x 2 — r 2 /3y 2 — r 2 
(3y 2 — r 2 /3x 2 — r 2 ) 00 pattern; dot-dashed line is evaluated 
within the linear chain. Inset: polarization for different values 
of U. 

WFs are bond centered, their centers located exactly in 
the middle of each bond connecting two neighboring sites 
along the zigzag chain, and ii) when 00 is established, 
WF centers move (in the same direction along the zigzag 
chain) toward lattice sites. This unveils the reduced sym- 
metry of the system, that loses the center of symmetry 
in the middle of each bond, thus allowing for a finite P. 
The displacement direction is related to the character of 
the occupied orbitals, in the sense that each WF center 
moves from the middle of the bond toward the neighbor- 
ing site where the occupied orbital is aligned parallel to 
the vector connecting each pair of sites. Since the 00 
pattern is 3x 2 — r 2 /3y 2 — r 2 (or 3y 2 — r 2 /3x 2 — r 2 ) and the 
direction along the chain changes as {..x,y,x,y..}, there 
is a net " leftward" (or " rightward" ) displacement of WF 
centers. 



Because of the analogies shared with the AFM-CE 
phase relevant for half-doped manganites, showing dou- 
ble zigzag spin chains [lU, let us briefly discuss the out- 
come of our analysis in that case. The Hamiltonian 
([2]) is still suitable to describe the motion of e g elec- 
trons within chains where the hopping direction changes 
as {..£, x, y, y.-}-, as expected on the background of the 
ferromagnetically aligned core-spin chains characteristic 
of the AFM-CE configuration. We can then distinguish 
between corner Mn sites, where hopping changes direc- 
tion, and bridge Mn sites, where it does not[2(|. As 
pointed out before, the system is a band insulator even 
in the absence of any interaction, due to the effective 
dimerization induced by the phase change in the hop- 
ping amplitudes at corner sites. However, an 00 pat- 
tern of 3a; 2 — r 2 /3y 2 — r 2 states already occurs on bridge 
sites [ljl |2(J, at odds with the AFM-E phase considered 
so far, where only corner sites appear. Turning on the 
local interaction gives rise to a charge transfer from cor- 
ner to bridge sites, but leaves the 00 pattern unchanged. 
We can evaluate P along the same lines described before. 
We notice however that electrons pick up a phase on cor- 
ner sites, depending on the change of hopping direction, 
and a different one on bridge sites, related to the orbital- 
ordered state, which do not interfere with each other. As 
a result, we find that P is always zero with a CE-type 
constraint on the electron motion. 
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